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Abstract. We generalize the Raychaudhuri equation for the evolution of a self gravitat- 
ing fluid to include an Abelian and non-Abelian hybrid magneto fluid at a finite temper- 
ature. The aim is to utilize this equation for investigating the dynamics of astrophysical 
C^) high temperature Abelian and non-Abelian plasmas. 
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1. Introduction 

I 

It is perhaps one of the finest tributes to the simplicity and elegance of Einstein's 
theory of Gravitation and to the human spirit, that a college teacher sitting in 
Kolkata was able to conceive the equation for the evolution of a gravitating fluid 
now known as Raychaudhuri's equation [1]. The equation served as a lemma for 
the Pcnrose-Hawking singularity theorems and for the study of exact solutions in 
general relativity [2,3]. It provides a simple validation of our expectation that grav- 
itation should be a universal attractive force between any two particles in general 
relativity [4], This equation has stood the test of time and has been generalized 
in many ways. It has found applications in modern theories of strings and mem- 
branes [5,6]. We attempt a more modest generalization to tackle the statistical 
properties of a hot astrophysical plasma. These may be electromagnetic or chro- 
momagnctic (quark gluon) plasmas. We are guided here by a recent formalism 
that has been used to investigate the dynamics of a hot charged fluid in terms of a 
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hybrid magneto-fluid and the changes brought about to the Raychaudhuri equation 
by the introduction of a statistical attributes associated with finite temperature [7] 
are many and interesting. We give also give a brief outline of how to generalize this 
to the evolution of gravitating non-Abelian plasmas in the early universe [8]. 

2. Raychaudhuri's equation for charged gravitating fluids. 

In this section, for the sake of completeness, we first review the derivation of the 
standard Raychaudhuri equation for the evolution of a self gravitating fluid. We 
follow the derivation from Raychaudhuri's original work for geodesic flow [9]. We 
then incorporate the Lorentz force into the Raychaudhuri equation for non geodesic 
flow in the presence of an electromagnetic field. 

Consider a self gravitating fluid in which the fluid particles have a velocity U^. In 
a spacetime with metric g^, these velocity vectors allow us to define an orthogonal 
projection vector 

= g^ + UftUv. (1) 

The following properties hold for h^ v : 

• 1. W — g^ + U^U^ — 3, assuming V^V^ = — 1. This implies that is three 
dimensional. 

• 2. 

tf"V = U *9^ + U^U V = (2) 

and 

K V U V = 9lxv U v + U^U V U V = U^-U^ = (3) 

This implies is orthogonal to h^. Since is tangent to timelike curves, 
hfj, v is purely spacelike. It is the metric of the three dimensional hypersurface 
orthogonal to U^. 

For any time-like vector in a spacetime, the geodesic deviation equation can 
be written as 

U^ap - U".,f, a - W aa pU° (4) 

Although Eqn.4 is the geodesic deviation equation, it holds for any timelikc 
vector and not just the tangent vector to the geodesic. This is an important 
observation since Electromagnetic and Yang-Mills fields cause the fluid motion to 
be non-geodesic. Let us consider the parallel transport of B a p — U a[ p. It has 
the following properties: 

. a.)U a B aP = U a U a , fi = \(U a U a ), fj = 0. 

• b.)B a f3ljP = Ua-pljP 7^ unless geodesic. 
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Thus, B a p is a purely spatial tensor for geodesic motion, but fails to be fully spatial 
for non geodesic motion. In particular, with an EM field present 

U a .,pUP = ^-F a0 U^ (5) 
m 

Following Raychaudhuri, let us decompose B a p into its irreducible parts: 

• i. The antisymmetric part ; f2 a /3, 

• ii. The symmetric, trace free part ; S a /3, 

• iii. The trace ; \ < dh a0 . 
We can write 

B a f} = -Ua^W'Up + S Q /3 + Qa/3 + \®h a {}, (6) 

with 6, giving the rate of dilation of a three space element locally orthogonal to 
U a , emerging as the trace 

B a a = U a , a = 6. (7) 

We can define 

o Q/3 = htu a .p - u (i , a ) + \{u a u t} - upu a ), (8) 

where, U a = U a .pU^ is the acceleration vector, which gives the departure of the 
velocity field from geodesity. Similarly 

s Q/3 = ^{u a ,f) + u fi , a ) + \{u a u + UpU a ) - leh aP . (9) 

Consider 

u a n al3 = ^{u a u a , (i - u a u p , a ) + ^(u a u a u fj - u a u p u a ) 

= -\{u a ., (i u a + \{u , a ) + ^(u a u a . n ir")Up + l -u fi , a u a = 0, (10) 

and 

n^u - \{u a , p uP - u p , a u p ) + ^(li^Ufj - u p u a u p ) = o. (11) 

These relations imply that Q a p is spatially antisymmetric. Similarly, 

U a ^ a(i = \{U a U a .jj + U a U (j , a ) + ^(u a u a ii fj + U a U fi U a ) - ^W a h afi 

= \{Uf>-Up) = 0. (12) 
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and, 

^ aP u t3 = \{u a .,pU t3 + u fj , a u fj ) + \{u a u fi u fj + u fi u f 'u a ) - \@KpU^ 

= \{U«-U a )=0. (13) 

So we have taken care of the non geodesic nature of U a to define T, a p and Cl a p, by 
including the acceleration due to the EM force and made them spatial. 

Let us now consider a parallely transported B a p, B a p ;l = U a; pj. Eq. 4 allows 
us to write 

-Bq/3; 7 = (C^a; 7 /3 - Ran-fpU^) (14) 

Thus, 

= (U a , 7 U^) ; p - U a „lP ;P - Rv^pWlP 
= {U a );p - B^Wp - R a ^pU»U">. (15) 
Taking the trace over a and (3 we get 

B a ar( U~< = ([/<% - B ai B^ - R a ^ a U»U\ (16) 

But 

s a a = e (17) 

since Yi a p is traceless and Q a p is antisymmetric. Therefore Eq.16 becomes 

9 ;7 t/ 7 = (U a ). a - B a , ( B~< a - R^U^IP (18) 

where i? M7 is the Ricci tensor. Using the Lorentz force law equation U a — ^F a 7 U" r , 
we have 

= (— F"C/ 7 ). a - B a ^ a - R^U^U" 
m 

= + —F a ~B^ a - B a . t B-< a - R^U^U" 

m m 

= -F a rQ C/ 7 + —F a ^W a - B a ^B~< a - R^U^XF (19) 
m m 

Where we have the antisymmetric nature of F a p, so that F a pB^ a — F a pfl al3 . We 
also have 



B al B~« a = S Q7 S 7Q + n aj W a + -6 2 /i a7 /i 7a (20) 



1 
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The Ricci tensor is related to the energy momentum tensor by the Einstein's equa- 
tion 

Rv* - \R9^u = ( 22 ) 
or, equivalcntly, by 

R„ v = 87r(t M „ - l -t gilv ) (23) 
Substituting for t we find 

e = — F% a m + ^-F a ^w a - - o, a ^w a 

m m 
- ie 2 - 8tt(V - \t 9llv )U»U v . (24) 

With our notation, (U^U^ = — 1), the limit — > is the standard (geodesic) 
Raychaudhuri equation. T, a fj is called the shear tensor, £l a p is called the vorticity 
tensor and 9 is called the expansion scalar. These names arise because one can define 
an orthogonal triad of vectors in the three dimensional hypersurface orthogonal to 
U^. Under Fermi- Walker transport, the volume defined by the triad expands, gets 
deformed and also rotated. 



3. Raychaudhuri's equation for a "unified" charged gravitating fluid 

For the use of the Raychaudhuri equation to describe astrophysical plasmas in which 
gravitational effects are compatible with the fluid attributes, we have to account 
for the temperature and pressure of the fluid. For this, we incorporate into a small 
volume element of the fluid a statistical factor / which represents a temperature 
dependent statistical attribute of the fluid, and is related to the enthalpy h, the 
scalar density in the rest frame n and the mass m of the fluid particles by the 
relation h = mnf(T). When one does the kinetic theory of high temperature 
plasmas f(T) seems to emerge as the most useful variable to represent temperature 
effects. For rclativistic plasmas h — ™t ^(m"/?) an< ^ /(-^) ^ s P urc ly a function 
of temperature. The velocity vector of the fluid is obtained as the average velocity 
of this small volume of the fluid and is written as = fU^. This drastically 
alters the character of the terms in the evolution equation. For example, now 
V% = -f in contrast with U^U^ = -1 and unlike U a U a . fj = ±(U a U a ) ;0 = 
0,, now V a V a[ p = —fdpf. These terms significantly change the spatial terms of 
the Raychaudhuri equation and necessitating a generalization to account for these 
statistical factors. Indeed, this statistical limit for the fluid velocity [7], unlike the 
particle limit, provides a natural factor for producing acceleration forces from within 
the fluid due to pressure and temperature gradients. In the unified magnetofluid 
picture one can write the equation of motion of a magneto-fluid with entropy er as 

T3 v g = gM^Un (25) 
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where 



771 

= + —S^, (26) 



and 

Spv = dp{fU v ) - d v (fU„), (27) 

are antisymmetric second rank "flow" tensors defined in [7]. In this sense, M^ v 
represents an anti-symmetric " unified" field -flow tensor constructed from the kine- 
matic (U^) , statistical (/(T)) and electromagnetic (F^) attributes of the magneto- 
fluid. 

While the statistical factors provide us with additional acceleration terms which 
alter the purely spatial character of the shear and vorticity tensors, there are some 
relations that remain unaltered. A trivial example is that V* remains orthogonal 
to the hypersurface with metric h a p defined as the projector h a p = g a p + U a Up. 
To derive the geodesic deviation for notice that 

V^afi = (/!7% + U^f, a ), 

= f-,pU», a + U^,pf, a + U»f ;a/3 + fU" iaP . (28) 

It is easy to see that anti symmetrization in indices a, (3 reduces the geodesic devi- 
ation to 

= WaapV, (29) 

and its character is unchanged. 

The question that arises now is if it is possible to define the generalizations of 
the standard definitions of the shear and vorticity tensors in a similar manner, and 
can they be constructed to be purely spatial tensors. 

Let us, following the standard procedure, decompose B a p into its irreducible 
parts and define: 

B a0 = S Q/3 + h af} + ^@h a p + jV a dpf. (30) 
Let us examine the trace of B a p as this defines the expansion scalar O 

B a a = e + jV a d a f = v a ia (31) 

where we have assumed that the generalized shear tensor T,^ is symmetric and 
traceless. If we now define 

S Q/3 = ^04 ;/3 + Vp, a ) - ^Gh af} - ^j(V a dpf + Vpd a f), (32) 
then its trace 
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s% = v ; „ - e - iva^/ (33) 

because of equation(31), goes to zero; the trace free condition also reproduces the 
required definition of 0. 

The generalized vorticity tensor Q a p may also be written as 

n a p = \{v a -,p - v . a ) - ^j(v a d p f - Vpd a f). (34) 

We see that the tensor S^ v defined in Eqn.27 allows us the following identification: 

= ^s 0a - ^j(v a d fJ f - Vf,d a f). (35) 

Then, a straight forward analysis shows that 

V a {t aP + n a(i )=0. (36) 

An easy consequence of this result is that 

V a B aSi = 0. (37) 

From earlier work, we know that the gradient of /, for a perfect fluid, is related to 
the pressure gradient and to Lorentz forces from the presence of a electromagnetic 
field [7]. So, unlike the standard particle picture of the shear and vorticity, the 
unified picture of shear and vorticity automatically includes accelerations coming 
from the internal forces of the fluids. This in turn ensures that the shear is now not 
a purely spatial tensor; the vorticity tensor behaves in a similar manner. Indeed, 
we see from the above equation that it is the sum of the shear and vorticity tensors 
that allows simplification. 

Let us examine BapV 13 : Again, it is easy to see that 

{t aP + + jV a dpf)V p = V a , V = V a , (38) 
implying B afj 

V a B afi = (39) 

and 

BapVP = V a .pV fj = V a . (40) 

The decomposition of V a - t p into its irreducible components can therefore be written 
as: 

V a .f} = B a(} = t a[j + tt a(} + ^Gh aP + jV a dpf. (41) 

So, while the generalized expansion scalar, shear and vorticity tensors show evidence 
of internal fluid forces, the tensor, B a p — V a -p itself shows similarity with the 
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original tensor, B a p — U a -p denned earlier in the particle picture of the gravitating 
fluid. I.e; 

V a B afi = 0, (42) 

and 

B aP VP = V a .,f,V0 = V a . (43) 

A careful analysis shows why this is so. From the definitions of 0, T,^ and (l^ u , it 
is easy to see that the following equations relate the scalar of expansion, shear and 
vorticity in the particle view to those in the unified view of gravitating fluids: 

e = /e 

*V = (44) 

which then imply, 

B„ v = fB^ + V^OJnf. (45) 

In analogy with the derivation of the standard Raychaudhuri equation, if we define 
B a f3 = V a ;p, then from the deviation equation (29), 

= {(v a .„v>). tf) - v a .„v>. tf) - R a a~ lf) v°v^) 

= (V a );p - B ai B\ - R aaiP V a V) (46) 

Taking the trace over the indices a/3 and using B a a = V a - a = + jV a d a f, the 
new expansion scalar, we have 

(0 + jVdnfinV = (V". a - B ai W a - R^VV). (47) 

As before, we use Einstein's equation to get 

(0 + jVdnfinV = (V-. a - B aP & 3a - 8tt(V - ^tg^VV). 

(48) 

The acceleration forces are twofold, one coming from the em forces within the 
fluid and the second from the fluid forces themselves (pressures, etc). To examine 
what these are, we return to the equation of motion for a fluid with entropy a in 
the unified picture Eqn.26 and use the fact that 

V - W(T)U V ) - d v (f(T)U v ) = - V KV (49) 

For an isentropic fluid, 
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U^d^a = o-pV* = 



and therefore, we can write 



Td v a = -TU»[U^d v a - U v d^d 
T 
1 



= —nV^dua-Vud^a] 



Putting this together with (26), we find 
or, 



—V^V^a - V v d^\ = qU»(F vll + -S„„) 



Therefore, 



a mT 



and, 



a mT 
V v = -fd v f + ^(V»F^ + —V^8 v a - V v d^\) 

Defining a pure fluid factor, as 
T 

= -j{<j-hV v - a^Vp), 
the acceleration vector simplifies to 

V v = -fd v f + -V»{F Vil + -N vv ). 

m q 



Defining 



r - F 4- — N 



substitution into (48) yields 

(6 + -V^d.f),^ = (-.fd a f + —VpG a(S ). a - B a f 3 B (3a 
J m 

- tort^VV - An ft 

which simplifies to 
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+ C = (-fd a f)-, a + ?-Vp ;a G a P + —VpG a ® - a 
m m 

- B a/3 B 0a - tort^VV - Airft (60) 

where we have defined V^d^lnf — (. Using the relationship between and B^, 
we find 

9 + C = (-fd a f)-. a + ^V fj .. a G ati + —Vf3G al3 - a 
m m 

- fB af ,fBP a - 2B a *V r d a lnf + ( 2 - tort^VV - Avft 

= —Vfs ;a G af3 + ±V p G a e. a - fB a pfB^ a 
m m 

+ ( 2 - f(d 2 lnf) - 8nt^V»V - 4nft (61) 

In the limit / — ► 1 as T — > this reduces to the standard Raychaudhuri 
equation . Simplifying the left hand side, 

6 = — V /3 . a G a0 + — VfjG a/3 . a - fB a pfBf> a + C 2 
m m 

-C-.f (d 2 lnf) - SttV W - ^ft (62) 
Simplifying the B^ V B U>1 terms, we have 

§ = ^-S a(3 G at3 + ^V (3 G ati , a - {t a pt? a 
2m m 

+ fi^fi** + ie 2 ) + C 2 - C - f(d 2 lnf) - 8wt^V^V v - Airft (63) 

This final form of the generalized Raychaudhuri equation shows all the terms 
which reduce to the standard equation in the limit f(T) — ► 1. 



4. A Raychaudhuri equation for the unified non-Abelian magneto fluid. 



In a recent paper [8], we had generalized the proposed unification of the Abelian 
magneto fluid [7] to a non-Abelian magneto fluid. The equations of motion for 
the magneto fluid clearly indicated the possibility of solitonic solutions which are 
normally absent in the Abelian case. The inherent non-linearities present in the 
Yang Mills magneto fluid allow such possibilities and a natural question arises 
as to what will self gravity do to the dynamics of a non-Abelian magneto fluid. 
This is particularly relevant in view of the compelling experimental evidence from 
the relativistic heavy ion collider at Brookhavcn National Laboratory (BNL) that 
the universe in its first few moments may have existed as a quark-gluon fluid. 
Since, large gravitational fields are also present in this epoch, it provides us with 
a motivation to give a generalization of the Raychaudhuri equation for the non- 
Abelian magneto fluid. 

The suggestion in [8] was that each worldline would now carry an internal index 
labelling the non-Abelian species. Generalizing the Lorentz force law to the non- 
Abelian case, we had derived the equation of motion for the magneto fluid in 
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terms of a unified antisymmetric tensor, M % ^ = F z ^ v + ^-S 1 ^, where F l ^ u is 
the standard Yang Mills field strength tensor while S 1 ^ is given by: 

S\ v = d^fU\) - d v {fU\) - igftAp, U v } 1 + igf[A v , - imf[U^, U v \ 

(64) 

with the gauge covariant derivative being defined by T>^ = — i<?L4 M ,]. The 
non-Abelian fluid equations of motion corresponding to Eqn.26 for a "Yang-Mills 
Magneto-fluid" , with entropy a are given by 

T3 v a = gM% v U atl . (65) 

For a non-Abelian magneto fluid velocity vector U 1 ^, the deviation equation can 
easily be written as 

XP*. va - U^, av = R% va U lp . (66) 

Examining = U l ^ u , 

From our earlier work [8], we write 

tr% = IP^Ui" = 9-F\ v \jr (68) 
m 

for the generalization of the Lorcntz force law in the particle picture of the 
non-Abelian magneto fluid; the trace is over the internal, gauge group indices 
(i = 1...N = dim(gauge group)). Clearly, since now, IP^U^ = —N, we shall 
henceforth, assume the U l M are normalized so that U l pU^ = — 1, i.e; we assume 
U\ — ► U\/VN. With this proviso, we can define an orthogonal projection tensor 
as before 

V = 9i»> + U\U iv . (69) 
It follows that 

U m K v = U^ 9llv + U^U\U lv = U j v - 5 ij U iv = (70) 

and 

h^U jv = g^U jv + U\U iv W v = XJ\ - S^U ifl = (71) 
where we have used 

U\U j » = -S ij (72) 
from which it follows that U % = —1. Equation(67) can be reduced to 

B jxv;aUi — ( — F jxcJJi ) ;l / — RppvaU P U { p U p-, a Ui - v . (73) 
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Tracing over indices /U and v, and defining 



e J - B\», (74) 

we find 

9 = {^F^ a un^ ~ R pa U ia U t p - B\ a B^ (75) 
Once again, we can decompose the tensor B z flv into its irreducible parts and write 

I, 

3 



B\ v = U\, v = Y,\ v + n\„ + -Q% v . (76) 



It is to be noted here that we are considering a special sector of the magneto fluid 
dynamics, that of a non- interacting sector of the full non-Abelian fluid as is seen 
from our definition of the acceleration vector a M = TJ % ^ v TJi v . From a fluid point of 
view, we are dealing with a multi-species model of the non-Abelian fluid. The full 
non-Abelian interactions of the fluid will be explored in a future study. Returning 
to the decomposition of B l ^ u , we write 



— -(TP - TP \ 
2L iii> 2 u '^' 



(77) 



Clearly, these definitions give a decomposition of B l ^ v into its irreducible compo- 
nents. The resulting Raychaudhuri equation, upon substitution of B l llv Bi Vfl is, as 
expected, a multi-species generalization of the standard Raychaudhuri equation. 
This Non-Abelian generalization of Raychaudhuri's equation gives us a means of 
studying the shear and vorticity of quark-gluon astrophysical plasmas. 



5. Conclusion 



High temperature plasmas have an important role to play in the early universe. 
In what is known as the " classical" or the " radiation" epoch of the universe, both 
gravitational and high temperature effects are of equal importance. Thus for the 
evolutionary dynamics of these hot gravitating plasmas, a generalization of the 
Raychaudhuri equation to include finite temperature fluid forces as well as elec- 
tromagnetic effects was called for. Using a unified magneto-fluid approach to con- 
struct a generalized Raychaudhuri equation, we have attempted to respond to this 
call. Non-Linear plasmas also find an application in cosmology through the (now 
compelling) evidence that a non Abclian, non linear quark gluon fluid existed in 
the early epochs of the universe. To deal with the evolution of this fluid when 
the gravitational effects are strong , we have laid the foundation of a generalized 
Raychaudhuri equation for the evolution of non Abelian gravitating plasmas. Most 
interesting plasmas involve collective effects which are non-linear even in the special 
relativistic situation and the inclusion of gravity can only lead to more intriguing 



12 



Pramana - J. Phys., Vol. xx, No. x, xxx xxxx 



highly non linear phenomena. The investigation of further implications of the ideas 
germinated in this paper is a promising avenue for future research. Amalgamation 
and unification of ideas that cut across barriers always enriches physics and we feel 
that any work of this kind is a fitting tribute to a man whose physics cut across 
international boundaries. 
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